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ABSTRACT 

We present a gravitational lensing and photometric study of the exceptional strong lensing 
system SDSS J1538+5817, identified by the Sloan Lens ACS survey. The lens is a luminous 
elliptical galaxy at redshift zi — 0.143. Using Hubble Space Telescope public images obtained 
with two different filters, the presence of two background sources lensed, respectively, into an 
Einstein ring and a double system is ascertained. Our new spectroscopic observations, performed 
at the Nordic Optical Telescope, reveal unequivocally that the two sources are located at the 
same redshift z s = 0.531. We investigate the total (luminous and dark) mass distribution of the 
lens between 1 and 4 kpc from the galaxy center by means of parametric and non-parametric 
lensing codes that describe the multiple images as point-like objects. Bootstrapping and Bayesian 
analyses are performed to determine the uncertainties on the quantities relevant to the lens mass 
characterization. Several disparate lensing models provide results that are consistent, given the 
errors, with those obtained from the best-fit model of the lens mass distribution in terms of a 
singular power law ellipsoid model. In particular, the lensing models agree on: (1) reproducing 
accurately the observed positions of the images; (2) predicting a nearly axisymmetric total mass 
distribution, centered and oriented as the light distribution; (3) measuring a value of 8.11 jloisg x 
10 10 Mq for the total mass projected within the Einstein radius of 2.5 kpc; (4) estimating a total 
mass density profile slightly steeper than an isothermal one [p(r) oc r~ 2 ' 33 -° - 20 ]. A fit of the Sloan 
Digital Sky Survey multicolor photometry with composite stellar population models provides a 
value of 20j^4 x 10 10 Mq for the total mass of the galaxy in the form of stars and of 0.9+q'3, for the 
fraction of projected luminous over total mass enclosed inside the Einstein radius. By combining 
lensing (total) and photometric (luminous) mass measurements, we differentiate the lens mass 
content in terms of luminous and dark matter components. This two-component modeling, 
which is viable only in extraordinary systems like SDSS J1538+5817, leads to a description of 
the global properties of the galaxy dark matter halo. Extending these results to a larger number 
of lens galaxies would improve considerably our understanding of galaxy formation and evolution 
processes in the ACDM scenario. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: individual (SDSS J1538+5817) — 
galaxies: structure — dark matter — gravitational lensing 
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1. Introduction 

Early-type galaxies host the majority of the 
baryo nic mass observed in g alaxies in the Universe 
(e.g., iFukugita et al. 1998t iRenzini 20061 ); hence, 
deciphering the processes that lead to their for- 
mation and the mechanisms that rule their subse- 
quent evolution is a key cosmological issue. For in- 
stance, it is still debated whether early-type galax- 
ies form at relatively high redshift (z < 2) as a re- 
sult of a global starbur st and then passive l y evolve 
to the present (e.g., Eggen et al. 19621: Larson 



19741 : lArimoto fc Yoshiilll987l : lBressan et all 1991 ) 
or whether they assemble from mu tual disruption 
of disks in merging events (e.g., iToomrd 11977 : 
White fc Reesl Il978l) . Information with which 
to distinguish these scenarios lies in the char- 
acteristics of galaxy dark-matter halos. How- 
ever, the lack of suitable and easily interpreted 
kinematical tracers, such as HI in spirals, has 
made comprehensive studies on the dark mat- 
ter compon ent in early-type g alaxies rather dif- 



ficult (e.. 



Thomas et al 



Bertin et al.l 1992t iSaglia et al 
20071 120091 



1992; 



In the last few years, strong gravitational lens- 
ing has allowed astrophysicists to make great 
progress in the understanding of the internal 
structure of early- type galaxies. Through lens- 
ing, it has become possible to address in detail 
some fundamental problems related to the mecha- 
nisms of formation of early- type galaxies, like the 
determination of the amount and distribution of 
dark matter (e.g., Gavazzi et al.ll2007 ; Grillo et al 



2008cL l2009t iBarnabe et al.l l2009h or the investi 



gation of the total ma ss density profile and its 



redshift evolution (e.g ., Treu fc Koopmansl 2004 : 



Koopmans et al.l 12006). Several algorithms have 



been developed in order to fit the observational 
data of a strong gravitational lens system and, 
thus, to reconstruct the properties of a lens. 

Simplifying, a first difference between codes is 
the use of a pa rametric model ( e.g., qravle nj^ 
Keetonll2001allbl: ISeitz et al.lll998l: IWarren fc Dve 



20031 : lHalkola et al.l 120061 : IRzepecki et al.l l2007t 
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Grillo et all l2008ch or a non-parametric model 
(e.g.. PixeLenM ISaha fc Williamsl20o"l iKoopmansl 
200,4 ISuvu et al.ll2009h to describe the mass dis- 
tribution of a lens. In the former case, the mass 
distribution of a lens is assumed to be accurately 
described by an analytical expression: the funda- 
mental scales of the model are determined by com- 
paring the observed and model-predicted proper- 
ties of the multiple images. In the latter case, 
a pixelated map or a multipole decomposition 
of the surface mass density of the lens is usu- 
ally estimated through a statistical analysis that 
requires, in addition to the observational informa- 
tion, some extra physically plausible constraints, 
called priors, on the surface mass density distri- 
bution (e.g., positive-definite and smooth) of the 
lens. On the one hand, parametric models provide 
a great deal of freedom and complexity, but they 
do not cover "naturally" all the possible realis- 
tic mass distributions (for instance, surface mass 
density distributions with twisting isodensity con- 
tours); on the other hand, even if non-parametric 
models are more general, their number of degrees 
of freedom is often much larger than the con- 
straints and this can result in three-dimensional 
density distributions that are dynamically unre- 
alistic or unstable. A viable solution to obtain 
physically significant density distributions is to 
consider a framework where the mass distribu- 
tion of the lens is reconstructed by combining 
in a fully self-consistent way both gravitational 
lensing and stellar dynamic s measurements (e.g., 
Barnabe fc Koopmansl 2007 ). 



Lensing codes are further distinguished by 
the fact that the multiple images and the cor- 
responding sources can be m odeled as point- 
like (e.g., qravlens ; PixeLens; Seitz et al] 19981: 
l Halkola et aTT 2006) or ext e nded (e.g.. I Warren fc Dvel 
20031: IKoopmansl [2005 IRzepecki et al.l 120071: 



brillo et al.ll2008cl ISuvu et al.ll2009h objects. In 
the context of point-like algorithms, the best- 
fit model is defined as that which minimizes the 
chi-square between the measured positions of the 
centroids of the images and the positions repro- 
duced by the model, weighted by the measure- 
ment uncertainties. Additional chi-square terms 
that quantify the agreement between the observed 
and model-predicted relative fluxes and time de- 



2 http: / / www.qgd.uzh.ch / projects / pixelens / 
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Fig. 2. — SDSS spectrum obtained within an aper- 
ture of 3" diameter centered on the lens galaxy. 



lays of the multiple images can also be included. 
For extended algorithms, the goodness of a model 
is estimated by comparing on a pixelated grid 
the image surface brightness morphology and dis- 
tribution which are observed to those which are 
reproduced by the model (after convolution with 
the relevant point spread function). 

The relative positions of a multiply imaged sys- 
tem can sometimes be measur ed with an accuracy 
of a few milli-arcseconds (e.g. JPatnaik fc Kemballl 
20011 ) and these positions represent the most im- 
portant constraints on the mass distribution of the 
lens. In fact, although the flux ratios of the mul- 
tiple images can be easily estimated and offer an- 
other important source of information, the sen- 
sitivity of the flux measurements to details such 
as the dark matter substructure of the lens, the 
extinction in the interstellar medium of the lens, 
the microlensing effects of the stars present in the 
lens, and the time variability in the source de- 
crease their potential. Time delays can also help 
to determine the mass distribution of a lens, but 
a statistically significant number of measurements 
of this kind is just starting to become available. 

The projected total mass enclosed within a 
cylinder of radius equal to the Einstein radius 
of a lensing system can be accurately measured 
by only fitting the observed position s of the mul- 
tiple images (e.g., Kochanek 199lL iGrillo et~al 



and rings associated with an extended source is 
necessary if the interest i s in the properties of 
both lens and source (e.g., ISwinbank et al. 2007 : 
Vegetti fc Koopmand 20091 ). By combining lens- 



ing and multiband photometric measurements, the 
amount of mass present in a lens galaxy in the 
forms of da rk and visible matter ca n be deter- 
mined (e.g.. IGrillo et al.ll2008al l2009h . 

In addition to projected total mass, lensing 
analyses allow one to estimate also th e total mass 



densi t y profile of lens galaxies (e.g., iRusin et al 



2003; IRusin fc Kochanek! 120051 ). This can be 



achieved either by combining in a statistical way 
lensing and stellar dynamics data in a sample 
of lens galaxies that are assumed to have a ho - 
mologous structure (e.g., Koopmans et al. 2006), 
by performing a joint lensing and extended stel- 
lar kinematics study in a single lens galaxy ( e.g., 
iBarnabe et aTll2Q09l : pTreu fc Koopmand 12004 ). or 
by using lensing only in exceptional lensing sys- 
tems that show multiple images of different sources 
probing wide angular an d radial ranges of the 
lens m ass distribution (e.g. JSvkes et al.lll998tlNair 



1998) 



In this paper, we study the lensing system SDSS 
J1538+5817, discovered by the Sloan Lens ACS 
(SLACS) survejU. This system is particularly in- 



2008c). whereas a detailed fit of the possible arcs 
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Fig. 3. — Color images of a 5"x 5" field around the gravitational lensing system SDSS J1538+5817, before 
(on the left) and after (on the right) the subtraction of an elliptical model, fitted on the luminosity profile of 
the lens galaxy. The images are obtained by combining the F606W HST/WFPC2 and the F814 W HST/ACS 
filters. 



Table 1: The lens galaxy. 
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References - iBolton et al.l l2008. 

Notes - Magnitudes are extinction-corrected modelMag (AB) from the SDSS. 



teresting because two different sources arc lensed, 
one into an Einstein ring with four luminosity 
peaks and the other into two images, by an early- 
type galaxy that has an almost circular projected 
light distribution. The large number of images at 
various angular distances from the galaxy center 
and the nearly perfect axisymmetric lensing con- 
figuration of the ring makes this system the ideal 
laboratory to disentangle the luminous and dark 
components of the lens mass distribution. 

The paper is organized as follows. In Sect. 2, 
we describe the observational data for the complex 
strong lensing system SDSS J1538+5817. We per- 
form parametric lensing analyses of this system in 
Sect. 3. Then, in Sect. 4, we investigate the lu- 
minous and dark matter composition of the lens 
galaxy. In Sect. 5, we summarize the results ob- 
tained in this study. Finally, in the Appendix, we 
model the lens mass distribution on a pixelated 



grid and compare these non-parametric results to 
those from Sect. 3. Throughout this work we as- 
sume H = 70 km s _1 Mpc -1 , fl m = 0.3, and 
Ha = 0.7. In this model, 1" corresponds to a 
linear size of 2.51 kpc at the lens plane. 

2. Observations 

The SLACS survey was started in 2003 and 
aims at studying, from a lensing and dynamics 
perspective, a statistically significant number of 
galaxies acting as strong lenses and located at red- 
shifts lower than 0.5. The candidate lenses were 
spectroscopically selected from the Sloan Digital 
Sky Survey (SDSSjl database by identifying those 
objects that show, in addition to the continuum 
and absorption lines of a possible lens galaxy at 



4 http:/ /www. sdss.org/ 
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Fig. 5.— The NOT/ALFOSC ID (on the top) and 2D (on the bottom) spectra of the lensing system SDSS 
J1538+5817. In the wavelength range shown here, the most prominent spectral features at the redshift of 
the lens (zi = 0.143) and sources (z s = 0.531) are marked, respectively, with dashed and dotted lines. The 
shaded regions indicate the principal telluric lines. The flux is given in arbitrary units. 




Fig. 4. — Observational configuration for the spec- 
troscopic measurements obtained at the NOT tele- 
scope with the ALFOSC instrument: position of 
the 1" long-slit. 
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Offset (arcsec) 

Fig. 6. — Intensity (in arbitrary units) of the three 
Hg, Hp, and [OIII] A5007 emission lines along the 
slit width. The dashed lines mark the positions 
of three observed peaks that trace closely the geo- 
metrical configuration of the D\, Qi, and Q4 im- 
ages. 
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redshift zi, one or multiple emission lines of a 
hypothetical source at a higher redshift z s . The 
most promising candidates were then observed at 
least once with the Hubble Space Telescope (HST) 
Advanced Camera for Survey (ACS) to confirm 
the lens hypothes i s (for further information, see 
Bolton et al]l2006ll2008h . This procedure resulted 
in the sample of 63 "grade- A " strong gravi t ationa l 
lensing systems presented in iBolton et al.l (2008). 
SDSS J1538 + 5817 is one of the lens galaxies dis- 
covered by the SLACS survey. The photomet- 
ric and spectroscopic observations taken by the 
SDSS are shown in Figs. Q] and [2 As described 
above, the redshifts of the lens galaxy and a source 
(zi = 0.143 and z s = 0.531) were measured from 
different spectral features. 

By making use of the publicly available ob- 
servations in the F814W and F606W filters of 
the HST/ ACS and Wide Field Planetary Cam- 
era 2 (WFPC2) respectively, we model the lu- 
minosity distribution of the lens galaxy and sub- 
tract the best-fit model from the images. In de- 
tail, a model for the lens galaxy is constructed 
by using an iterative procedure: first, presum- 
able background source images are masked and 
isophotal contours of the lens galaxy are derived 
for surface brightness levels separated by 0.1 mag 
arcsec -2 . Then, all isophote contours (even if 
partially mas ked) are fitted by ellip s es fol lowing 
the method of iBender & Moellenhofl (jl987t ). This 
provides five parameters (center coordinates, ma- 
jor and minor axis, and position angle) per sur- 
face brightness level. The resulting table of these 
parameters and associated surface brightnesses is 
employed to calculate a smooth elliptical model 
for the lens galaxy. This model is subtracted from 
the image, leaving as residuals only the images 



Table 2: Astrometric and photometric measure- 
ments for the multiple images. 
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of the background source galaxy. The residuals 
are used to check and improve the masking of 
the source galaxy. This procedure is repeated 
twice until the final model for the lens galaxy 
is obtained. In addition, images taken in differ- 
ent wavelength bands provide color information 
for lens and source which is of additional help to 
identify and separate lens and source components. 
Color composite images of the strong lensing sys- 
tem and the residuals after the lens galaxy sub- 
traction are shown in Fig. [3] 

The excellent angular resolution of the HST al- 
lows us to identify two systems of multiple images 
(for labels, see Fig. [7]): a double (D\ and D 2 ) and 
a quad (Qi, Q2, Qz, and Qi). The images of the 
two systems have different colors, but the average 
distances from the galaxy center of D\ and D 2 
and Qi, Q 2j Q31 and Qi are consistent within the 
errors. This fact would imply approximately the 
same redshift for the two sources, if the lens to- 
tal mass were close to an isothermal distribution. 
In addition, the absence in the SDSS spectrum of 
evident emission lines at a possible third redshift 
supports the hypothesis that the two sources are 
at the same distance to the observer. 

Since the precise knowledge of the redshift of 
the two sources plays a crucial role in the deter- 
mination of the total mass distribution of the lens 
galaxy, we decided to perform additional spectro- 
scopic measurements to understand whether the 
emission lines observed in the SDSS spectrum are 
associated to one or both of the lensed sources. 
The data were obtained on June 25, 2009 as a 
Fast- Track Observing Program (P38-428) with the 
Andalucia Faint Object Spectrograph and Camera 
(ALFOSC) at the 2.5-m Nordic Optical Telescope 
(NOT) on La Palma (Spain). We positioned a 1"- 
wide long-slit centered in Qi and passing through 
Di and Q 4 , as shown in Fig. We used ALFOSC 
with the 8 grism, that covers a wavelength range 
between 5825 and 8350 A with a dispersion of 1.3 
A per pixel. In good atmospheric conditions (see- 
ing between 0.7 and 1") and in the same obser- 
vational configuration, we obtained six exposures 
of 24 minutes each, resulting in a total integration 
time of 2.4 hrs. 

In Fig. [51 we show the wavelength-calibrated 
ID and 2D spectra. We identify several prominent 
absorption lines at redshift 0.143 and at least six 
secure emission lines (Hs, G, H 7 , H^, [OIII] A4959, 



6 



and [OIII] A5007) at redshift 0.531. At the emis- 
sion line positions, the presence of two intensity 
peaks, below and above the continuum, is visible 
in the 2D spectrum. In Fig. [51 we plot the inten- 
sity of the Hs, Hp, and [OIII] A5007 emission lines 
as a function of spatial position along the slit's 
cross-section. This corresponds to a representa- 
tion of three sections of the 2D spectrum at the 
emission line abscissas. For these three emission 
lines, we distinguish three intensity peaks located, 
with respect to the continuum (lying between Qi 
and Qi), at angular positions consistent with those 
of the three images D±, Qi, and Q4. The mea- 
surement of the same emission lines at the same 
observed wavelengths proves in a conclusive way 
that the two sources D and Q are equally distant 
from the observer. We remark that the intensity 
values of the emission lines shown in Fig. [6] are 
differently contaminated by the lens galaxy flux. 

In Table [TJ we summarize the photometric and 
spectroscopic properties of the lens galaxy: the co- 
ordinates (RA, Dec, zi), the minor to major axis 
ratio (<7l) and its position angle (# 9L , degrees east 
of north), and the SDSS multiband magnitudes (u, 
g, r, i, z). In Tabled we report the coordinates 
of the multiple images (x, y, z s ) and the adopted 
position uncertainty on the first two coordinates 
{5x,y)i the relative flux of the double system com- 
ponents (/) and the respective error (5/), and the 
distance of the images to the galaxy center (d). 

3. Strong gravitational lensing 

We address parametric (Sect. I3.1[) point-like 
modeling of the strong gravitational lensing sys- 
tem. We focus mainly on projected total mass 
and total mass density profile measurements. A 
comparison with the results obtained from non- 
parametric models is provided in the Appendix. 

3.1. Parametric models 



Gravlens (|Keetonll2001af) is a publicly- available 
lensing software that, starting from the measured 
observables of a strong lensing system, recon- 
structs the properties of a lens in terms of an 
adopted model that is defined by some relevant pa- 
rameters. By using this code, we perform a para- 
metric analysis in which we describe the total mass 
distribution of the lens galaxy in terms of either an 
elliptical de Vaucouleurs model (deV), or a singu- 



lar isothermal ellipsoid (SIE) model, or a singular 
power law ellipsoid (PL) model (fo r further detail s 
on the model definitions, see e.g. iKeetonlliffllbh . 
Both a deV and an SIE model are characterized by 
five parameters: a length scale b (corresponding to 
the value of the Einstein angle #Ein in the circular 
limit), the two coordinates of the center (xi,yi), 
the minor to major axis ratio q, and its position 
angle q . For the deV model, we fix the value of 
the effective angle (8 e ) to that shown in Table[TJ A 
PL model is more general than an SIE model. In 
particular, the former requires as an additional pa- 
rameter the value of the exponent 7 of the three 
dimensional density distribution p(r) oc r~ 7 (an 
SIE model is retrieved by setting 7 equal to 2). 
The convergence n{x,y) of a PL model, defined 
as the surface mass density of the model divided 
by the critical surface mass densi ty of the studied 



lensin g system (for definitions, see lSchneider et al 



1992), depends on the previous parameters as fol- 
lows 

k(x, y) oc — . (1) 



Due to the normalization used in the code, 
gravlens provides values of a length scale b that 
are related to the values of b by a function /(■) of 
the axis ratio q: 



b = bf(q) 



(2) 



Varying the parameters of the two adopted 
mass models and the positions of the sources 
[{xD,yD);{xQ,yQ% we minimize a chi-square \ 2 
function. This function compares first only the ob- 
served (see Tabled]) and predicted positions of the 
multiple images [deV(nf), SIE (nf), and PL (nf) 
models] and then also the measured (see Table 
[2]) and reconstructed fluxes of the double system 
[deV (wf), SIE (wf), and PL (wf) models]. In the 
latter case, the flux of the source imaged twice 
is an additional free parameter of the models. In 
our lensing analysis, we decide to neglect the flux 
constraints relative the quad system because the 
presence of the Einstein ring prevents us from 
separating accurately the individual components. 
For the multiple images, we assume position un- 
certainties equal to the size of one pixel of ACS 
(0.05") and flux uncertainties as reported in Ta- 
ble [2] and determined by considering the different 
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Table 3: The best-fit (minimum x 2 ) parameters of the different models. 
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0.04 


0.866 


147.5 
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1.15 
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1.05 


-0.02 


0.02 
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147.6 




2.00 


3.77 
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PL (nf) 


0.82 


-0.02 


0.03 


0.800 


147.9 




2.33 


0.36 
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PL (wf) 


0.71 


-0.01 


0.02 


0.820 


148.0 




2.47 


0.99 


3 



Notes - The notation (wf) and (nf) indicates, respectively, if the flux measurements of the double system 
are included in the modeling or not. 



level of contamination on the flux estimates by the 
surface brightness of the lens galaxy. 

To estimate the statistical errors in the param- 
eters characterizing each model, we perform 2000 
X 2 minimizations on simulated data sets. These 
are obtained by extracting the image positions 
and fluxes from Gaussian distributions centered 
on the measured values and with standard devi- 
ations equal to the observational errors reported 
in Table [H In addition, starting from the sets of 
optimized parameter values, we estimate the total 
projected mass Mf°*(< Ri) enclosed within seven 
different circular apertures of radii Ri. The first 
three radii are chosen as the projected distances 
from the lens galaxy center of the inner image of 
the double system, the "average" Einstein circle, 
and the outer image of the double system. The 
remaining four radii are given by the values of the 
midpoints of the three segments defined by the 
previous three points and a further point at the 
same distance from the outer double image as the 
first point is from the second one [i.e., (0.26, 0.52, 
0.75, 0.98, 1.28, 1.58, 1.84)"]. 

The best-fit (minimum chi-square) parameter 
and x 2 values of the different models are sum- 
marized in Table [3] For all the models, we 
find that the best-fit x 2 values are smaller than 
the corresponding number of degrees of freedom 
(d.o.f.). This implies that the reconstructed po- 
sitions of the images are angularly very close to 
the measured positions of Table [2] In this sec- 
tion, we concentrate on the results given by the 
one-component PL and, as a comparison with the 
results of previous studies, simpler SIE models and 
only in the next section we will address the two- 



component mass decomposition. 

In Fig. for the best-fit SIE (nf) model we 
show the reconstructed positions of the sources 
and the caustics, the observed and reconstructed 
positions of the images and the critical curves, 
and the Fermat po tential (for definition, see 



Schneider et al.lll992T ) with its stationary points. 



The inclusion of the fluxes of the double system 
does not change significantly the best-fit parame- 
ters of the models. Moreover, the projected total 
mass estimates, which are presented below, are 
not particularly sensitive to the flux constraints. 
For these reasons, in the following we will mainly 
concentrate on the properties of the models that 
omit the additional source of information coming 
from the fluxes of the double system, i.e., the SIE 
(nf) and PL (nf ) models. We note that the best-fit 
b values are on the order of 1", the typical distance 
of an image of the quad system from the center of 
the lens (see Tables [2] and [3]). The best- fit values 
of the lens center and ellipticity show that the 
lens mass distribution is well centered and aligned 
with the galaxy light distribution (see Tables [1] 
and [3]). In particular, we remark that the total 
surface mass of the lens is well approximated by 
an axisymmetric distribution. The best-fit values 
of the parameter 7 of the PL models suggest that 
the lens total density profile is slightly steeper 
than an isothermal one. We estimate maximum 
time delays of approximately 30 and 3 days for 
the double and quad systems, respectively (see 
Table |4] and Fig. [7]). Finally, as far as the po- 
sitions of the sources are concerned, the model 
predicted angular distance of the two sources is 
between approximately 0.5 and 0.7"(see Table [5]), 



S 



corresponding, respectively, to 3.2 and 4.5 kpc at 
a redshift of 0.531. 

In Figs. [8]and[9l we plot, respectively, the joint 
probability distributions of the SIE (nf) and PL 
(nf) model parameters, with the 68% and 95% 
confidence regions and the 68% confidence inter- 
vals. These intervals are determined by excluding 
from the 2000 \ 2 minimizations the 320 smallest 
and the 320 largest values for each model param- 
eter. We have checked that the error estimates 
determined in this way are unbiased and equiva- 
lent to the uncertainties provided by a full Markov 
chain Monte Carlo analysis. The comparison Figs. 
[8] and [9] shows clearly that adding the exponent of 
the total mass distribution among the parameters 
increases their degeneracies, hence their error es- 
timates. The probability distribution of the posi- 
tion angle 9 q is bimodal, with a secondary peak 
located nearly 90° away from the primary one, 
found at approximately 150°. From the last col- 
umn of plots in Fig. [8l we see that the secondary 
peak is included only in the 95% CL regions, and, 
from the last panel of the same figure, we note 
that the low values of q are associated with al- 
most circular models (i.e., q ~ 1). The bimodal 
distribution of the values of the lens position an- 
gle can then be explained by looking at the source 
plane of Fig. If the axis ratio of a model is close 
to one, the shape of the tangential caustic is ap- 
proximately symmetric with respect to the center 
of the lens and the surface enclosed by this caus- 
tic is small. In the same limit, the radial caustic 
is well approximated by a circle centered on the 
lens center. From these considerations, it follows 
that the expected positions of the images are al- 
most invariant under a rotation of 90° of the lens 
mass distribution (supposing the positions of the 
sources are fixed). 

The degeneracies between b, q, and 7 are con- 
nected to their relations defined in Eqs. |T]) and 



Table 4: The model-predicted time delays for the 
best- fit model parameters given in Table [3J 



Model 












(days) 


(days) 


(days) 


(days) 


SIE (nf) 


24.8 


0.66 


0.70 


2.88 


PL (nf) 


33.5 


0.82 


0.87 


3.26 



((2|). In particular, the strong anti-correlation be- 
tween the value of the length scale and the steep- 
ness (see the fifth panel of Fig. [9|) is caused by 
the fact that the Einstein ring of a circular lens 
model defines a region on the image plane within 
which the average value of the convergence k is 
equal to one. In order for this equality to be ap- 
proximately valid inside the average circle defined 
by the positions of our quad system, from Eq. (fl]) 
and by holding the value of q fixed, it follows that 
a higher value of b requires a lower value of 7, and 
vice versa. 

The previous considerations on the almost 
model-independent average value of n inside the 
Einstein ring can also be translated in terms of to- 
tal mass estimates within the same ring. Distinct 
models, defined by different parameters, that can 
reproduce well an approximately complete Ein- 
stein ring, provide total mass measurements in- 
side this typical aperture that differ by only a few 
percent. This is shown in Fig. 1101 There, we 
plot the median values and the 68% confidence 
intervals (obtained by excluding from the 2000 
X 2 minimizations the 320 smallest and the 320 
largest mass estimates) of the lens projected total 
mass within the Einstein ring and measure values 
of 8.35t£?8 x 10 10 M Q for an SIE (nf) model and 
8.1lt f 9 x 10 10 M Q for a PL (nf) model. We notice 
that the median values of the 2000 Monte-Carlo 
cumulative total mass estimates do not necessar- 
ily follow a global PL model, but they have in 
principle more freedom. In fact, even if the total 
mass values of each of the 2000 models do follow a 
power law model precisely at all radii, the median 
values shown in Fig. [10] and used in the following 
for the luminous and dark matter decomposition 
are more general and do not provide the same 
value of the steepness 7 at each radial position. 
In general, for the two different models the total 
mass estimates, that are measured within various 
apertures (approximately between 1 and 4 kpc 



Table 5: The model-predicted source positions for 
the best-fit model parameters given in Table [3] 



Model 


{xd,Vd) 


(xq,Vq) 


d.D,Q 




(",") 


("/') 


(") 


SIE (nf) 


(0.27,0.46) 


(0.00,0.00) 


0.53 


PL (nf) 


(0.39,0.60) 


(0.02,-0.01) 


0.72 



9 



X 1.1 



i-xjoint probability distribution b-y, joint probability distribution b-q joint probability distribution b-S q joint probability distribution 



-0.15 -0.10 -0.05 -0.00 0.05 0.10 0.15 -0.15 -0.10 -0.05 -0.00 0.05 0.10 0.15 0.6 

x, (arcsec) y l (arcsec) 



60 80 100 120 140 160 
0, <<ieg) 



0.15 
0.10 
0.05 
-0.00 
-0.05 
-0.10 



x—y t joint probability distribution 



0.15 
0.10 
0.05 
-0.00 
-0.05 
-0.10 
-0.15 



x r1 joint probability distribution x— 9 joint probability distribution 

' ~ ; - 0.151 ' ' I 



0.10 
0.05 
-0.00 
-0.05 
-0.10 
-0.15 



-0.15 -0.10 -0.05 -0.00 0.05 0.10 0.15 0.6 

y, (arcsec) 



0.8 0.9 

1 



60 80 100 120 140 160 
0, <<ieg) 



0.15 
0.10 
0.05 

1 -° 00 

-0.05 

-o.io 

-0.15 



yi~ q joint probability distribution joint probability distribution 




0.15 
o.io - 

0.05 - 

o.oo j' 
0.05 - 

o.io - 

•0.15 L 



0.7 0.8 
1 



60 80 100 120 140 160 
O.OJeg) 

q-@ joint probability distribution 

i.0«7 




60 80 100 120 140 160 
0, (<!eg) 



Fig. 8. — Estimates of the errors and correlations in the parameters for an SIE (nf) model. Results of the 
X 2 minimizations on 2000 Monte-Carlo simulated data sets. Thick bars on the co-ordinate axes and contour 
levels on the planes represent, respectively, the 68% confidence intervals and the 68% and 95% confidence 
regions. For each model parameter, the 68% confidence interval is determined by excluding from the 2000 
X 2 minimizations the 320 smallest and the 320 largest values. 



from the lens center), are consistent, given the 
errors. We remark that fixing the exponent of the 



total mass profile (i.e., 7 equal to 2 for the SIE 
models) result in significant smaller uncertainties 
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Fig. 9. — Estimates of the errors and correlations in the parameters for a PL (nf) model. Results of the % 2 
minimizations on 2000 Monte-Carlo simulated data sets. Thick bars on the co-ordinate axes and contour 
levels on the planes represent, respectively, the 68% confidence intervals and the 68% and 95% confidence 
regions. For each model parameter, the 68% confidence interval is determined by excluding from the 2000 
X 2 minimizations the 320 smallest and the 320 largest values. 



in the total mass values. As mentioned above, by 
modeling also the fluxes of the double system we 
find total mass measurements that are consistent 



within la with the estimates obtained by fitting 
the image positions only. 
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Fig. 7.— Best- fit SIE (nf) model. Top left: source 
plane with caustics. The predicted source posi- 
tions of the double and quad systems are repre- 
sented by a plus and a cross symbol, respectively. 
Top right: image plane with critical curves. The 
observed and predicted image positions of the dou- 
ble (diamond and plus symbols, respectively) and 
quad (square and cross symbols, respectively) sys- 
tems are shown. Bottom left: contour levels of the 
Fermat potential for the double system. The im- 
ages are one minimum (Di) and one saddle point 
(D2). Bottom right: contour levels of the Fermat 
potential for the quad system. The images are 
two minima (Qi and Q3) and two saddle points 
(Q2 and Q 4 ). 



We generalize our result by emphasizing that 
the adoption of an isothermal model for strong 
lenses often provides a good fit of the observed 
images, but the errors on the projected mass es- 
timates may be considerably underestimated al- 
ready at projected distances from the center of the 
lens that differ from the Einstein radius by half its 
value. This fact has non-negligible consequences 
on the inferred properties of a lens dark matter 
distribution (see below). 

We notice that the value of 205 ± 13 km s^ 1 
for the central stellar velocity dispersion do , which 
is determined by rescaling the value of the SDSS 
spectroscopic stellar velocity dispersion measured 
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Fig. 10. — Comparison of the projected total mass 
estimates for the SIE (nf), SIE (wf), and PL (nf) 
models. For each aperture, the la confidence 
intervals are determined from 2000 Monte-Carlo 
simulations by excluding the 320 smallest and the 
320 largest mass estimates. 
The arrows show the projected distances of the 
observed multiple images from the lens center. 

inside an aperture of 1.5" [<r = (189 ± 12) km s _1 ] 
to an aperture of radius equal to 9 e /8, is consis- 
tent, within the errors, with the value of 215 ± 5, 
which is obtained by converting the total mass es- 
timates shown in Fig. \W\ for the SIE (nf) model 
into an effective velocity dispersion ctsie- 

We remark that the best-fit parameters of our 
SIE point-like models are consistent, given the er- 
rors, with the best-fit para meters of the S IE ex - 
tended model measured by iBolton et al.l ( 2008) . 
We also n o te that previous studies dKochanekl 
19931 Il994t iTreu et alllioOrl lOrillo et alJ l2008bh 
agree on finding that the central stellar velocity 
dispersion of early-type galaxies is a good estima- 
tor of the velocity dispersion of a one-component 
isothermal model. 

4. Luminous and dark matter 

We combine the surface brightness distribution 
measurement obtained from the HST images (see 
Sect. 2) with the multicolor photometric obser- 
vations of the SDSS (see Table 1) and the pro- 
jected total mass estimates determined from the 
lens modeling (see Sect. 3) to study the amount 
and distribution of luminous and dark matter in 
the lens galaxy. 

First, we compare in Fig. [Tl]thc surface bright- 
ness and the total surface mass [for the SIE (nf) 
model] isodensity contours of the best-fit models 
described in the previous section. We use two im- 
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Fig. 11. — Light and mass distributions. Isoden- 
sity contours of the best-fit surface brightness (on 
the left) and parametric [SIE (nf), on the right] to- 
tal surface mass profiles. The observed image po- 
sitions of the double (diamond) and quad (square) 
systems are shown. 

ages with the same area and pixel size, normalize 
the images to the sum of the values of all their 
pixels, and plot the same contour levels in both 
images. We observe that the distributions of light 
and total mass from the best-fit SIE (nf) model 
are nearly axisymmctric, but the former is slightly 
more concentrated than the latter. This can be 
inferred by looking at the positions of the inner 
and outer contour levels. The position angle of 
the surface brightness and total surface mass dis- 
tributions are consistent within the errors. Thus, 
we conclude that the light distribution is approx- 
imately circular symmetric in projection and it is 
a good tracer of the total mass distribution. 

Next, we fit the lens spectral energy distribu- 
tion (SED), consisting of the SDSS ugriz mag- 
nitudes (see Table 1), with a three-parameter 
Bruzual & Chariot composite stellar population 
(CSP) model computed by adopting a Salpeter 
initial mass function (IMF) and solar metallic- 
ity (for further details, see Grillo et al. 2009). 
The best-fit model, shown in Fig. [T^l provides a 
photometric (luminous) mass M*-^^ of the lens 

of 20±4 x 10 10 M Q . We then estimate the value 
of the mass in the form of stars •^*j 10 ^(^ R)> 
at a projected distance R from the center of the 
lens, by multiplying M*-^ Q ^ by an aperture fac- 
tor iap(< R), that represents the fraction of light 
measured within a circular aperture of radius R 
divided by the total light of the galaxy. The quan- 




T = 11.50 Gyr, t = 1.8 Gyr, W ^ = 2.0e+11 M 
153812.92 + 581709.8 , 2, = 0.143 
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Fig. 12. — SED and best-fit model of the lens 
galaxy SDSS J1538+5817. The observed total flux 
densities, measured in the u, g, r, i, and z pass- 
bands, and their la errors are represented by cir- 
cles and error bars. The best-fit is obtained by 
using Bruzual & Chariot 2003 models. On the 
bottom, the best-fit values of the age (T), the char- 
acteristic time of the SFH (r), and the luminous 
mass (M* hot ) are shown. 



tities introduced above are explicitly defined as 



phot 



(< «)=^phot/ap(<i*) 



and 



/ap(< R) 



J R I(R)RdR 



J™I(R)RdR 
where I(R) is the de Vaucouleurs profile 



I(R) = Iq exp 



- 7.67 



(3) 
(4) 

(5) 



with R e equal to D Q \ 8 e . As discussed above, 
the circular symmetry of the light distribution as- 
sumed in the previous three equations is a plausi- 
ble approximation for the lens surface brightness. 
In Fig. 021 we plot at different radii the projected 
total and luminous mass estimates obtained from 
the best-fit SIE (nf) and PL (nf) models of Sect. 
3 and the best-fit SED model. 

In the same figure, we show the fraction of pro- 
jected mass in the form of stars 



/*(< R) 



M; hot (< R) 
<!(< R) 



(6) 
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Fig. 13. — Comparison of the projected total mass 
estimates of the SIE (nf) (on the top left) and PL 
(nf) (on the top right) models with the luminous 
mass measurements from the best-fit SED model. 
The fraction of projected mass in the form of stars 
(on the bottom left) and the total and stellar mass- 
to-light ratios (on the bottom right) are obtained 
by using lensing and photometric information. In 
each panel, the arrows show the projected dis- 
tances of the observed multiple images from the 
lens center and the curves the la confidence inter- 
vals. 



and the total mass-to-light ratio 



M 



tot 



<<R) 



L B {< R) 



(7) 



plotted versus the projected radius R. We com- 
pare this last quantity with the values of the stel- 
lar mass-to-light ratio estimated from the best- 
fit SED model (M* hot L~ l = 5.0+U M Q L^ B ) 
and the evolution of the Fundamental Plane 
[M^Lg 1 = (6.1 ± 1.8) M Q L~^ B ] (for more in- 
formation, see Grillo et al. 2009). 

The need for a dark component to be added to 
the luminous one to reproduce the total mass mea- 
surements of the SIE (nf ) models is suggested by 
looking at the outer galaxy regions probed by lens- 
ing. Due to the larger error bars, the evidence on 
the presence of dark matter is reduced if the total 
mass estimates obtained from the PL (nf) models 



are considered. According to all the lensing mod- 
els (i.e., deV, SIE, and PL), a value of 0.9+^ for 
the fraction of projected mass in the form of stars 
over total is estimated at a projected distance from 
the galaxy center of approximately 2.5 kpc, and at 
4 kpc from the galaxy center a value of one for the 
same quantity is excluded by the SIE (nf) mod- 
els at more than 3a level. Moreover, at the same 
distance, the value of the total mass-to-light ratio 
determined from the SIE (nf) mass measurements 
is not consistent with the value of the mass-to- 
light ratio of the luminous component estimated 
from the galaxy SED modeling. Between 1 and 4 
kpc, the same decrease of /*(< R) and deviation 
of M tot Lg 1 from Lg 1 are also indicated 

by the values of the PL (nf) mass estimates, but 
these results are not highly significant because of 
the large uncertainties. 

Finally, by taking advantage of the total mass 
measurements available at different distances from 
the center of the lens (not only in the vicinity of 
the Einstein angle, as in the majority of the known 
lensing systems) , we decide to investigate the dark 
matter component in greater detail. To make pos- 
sible a direct comparison of our results with those 
obtained from dynamical analyses or cosmological 
simulations, we consider two-component models in 
which the luminous ph(r) and dark Pd(?") matter 
density distributions are parametrized by 

(3 - lL )M L r L 



Ph{r) = 
Pn(r) = 



47n-7L(r + r L ) 4 "TL 
(3-7 D )M D r D 



(8) 



47r r 7D( r + r D ) 4 -fD ' 

where Mj/ D is the total mass, r L / D a break ra- 
dius, and 7l/d t ne inner density slope of the lumi- 
nous and dark matter distributions. The density 
profiles of Eq. ||5J) are projected along the line-of- 
sight to give the corresponding surface mass den- 
sity profiles E L / D (i?): 



El(R) = 2 / 

JR 

£dCR) = 



Ph{r)r dr 



V 'r 2 - 


- R 2 


Pr>(r 


r dr 



R 2 



(9) 



which, once integrated, result in the following cu- 
mulative mass distributions A/l/d(< R)' 

M L {<R) = [ S L (i?)2vri?di? 

j0 
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M D (< R) 



E u (R)2irRdR. (10) 



The total density pr{r), surface mass density 
Et(-R), and cumulative mass Mr(-R) distributions 
are denned as the sum of the luminous and dark 
contributions 

Pr(r) = Ph(r)+pv(r), 
E T (i?) - £ L (i?) + £ D (i?), 
M T (< R) = M L (< R) + M D (< R) . (11) 

We notice again that the circular approximation 
is plausible for this particular lens. 

The luminous quantities introduced in the 
above equations are completely determined from 
the photometric observations. In fact, for the lu- 
minous component we have estimated the total 
mass Ml by modeling the SED and, to obtain a 
surface brightness profile close to a de Vaucoleurs 
profile, we assume a Hernquist (1990; 7l = 1 and 
r L = H e /1.8153) or a Jaffe (1983; 7 L = 2 and 
tl = i? e /0.7447) density profile. Then, we con- 
struct a grid of 13671 models for the dark com- 
ponent. The total mass Md, the break radius td, 
and the inner density slope 7d can assume values 
included between 0.1 and 100 times Ml, 0.1 and 
10", and 0.5 and 2.5, respectively. The first two 
intervals are divided logarithmically into 31 and 
21 points respectively, the last one linearly into 21 
points. The best-fit dark matter profile is found 
by minimizing the following chi-square function: 



X 2 (M D ,r D ,7 D ) 



E 



M T (< Ri) 



(12) 

In order to estimate the errors in the best-fit pa- 
rameters, we perform 500 Monte-Carlo simula- 
tions varying the total mass of the luminous com- 
ponent according to the corresponding measure- 
ment errors and the luminous break radius by as- 
suming a realistic 10% uncertainty. 

In Fig. 1141 we show the luminous and dark mass 
decomposition obtained from the best-fit (mini- 
mum chi-square) model and in Fig. [T5] the pa- 
rameter joint probability distributions. We decide 
to plot the best-fit dark matter model obtained 
by assuming a Jaffe profile (no significative differ- 
ences are present if a Hernquist profile is adopted) 
for the luminous component and considering the 



projected total mass estimates coming from the 
PL (nf) models. The confidence levels on the pa- 
rameter space of the dark matter component are 
expressed in terms of the luminous mass fraction 
/l = Ml/ (Ml + Md), i.e., the mass in the form 
of stars to the total mass of the galaxy, the ratio 
of the dark to luminous break radius r^/r^, and 
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Fig. 14. — Best-fit (minimum chi-square) lu- 
minous and dark matter decomposition, deter- 
mined by assuming a Jaffe profile for the three- 
dimensional luminous density and projected total 
mass measurements as estimated from the PL (nf ) 
model. The arrows show the projected distances 
of the observed multiple images from the lens cen- 
ter. 

We find a best-fit x 2 value of 0.8 with two de- 
grees of freedom (derived from the total mass mea- 
surements at the five central radii fitted by three- 
parametric models). We measure that the values 
of the dark matter density overcome those of the 
luminous matter density at radii larger than ap- 
proximately 1.5 times the effective radius of the 
galaxy (R e — 4.0 kpc). As in the previous sections, 
a three-dimensional total density profile close but 
not exactly equal to a function decreasing as 1 /r 2 
(i.e., an isothermal profile) is found. We note that 
the uncertainties in the dark matter parameters 
determined by using the projected total mass esti- 
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Fig. 15. — Estimates of the errors and correla- 
tions in the parameters related to the dark matter 
component: the luminous mass fraction /l, the 
dark to luminous break radius ratio r^/r^, and 
the dark matter inner density slope 7d- The pro- 
jected total mass measurements of the PL (nf) (on 
the left) and SIE (nf) (on the right) models are 
used. The small squares on the three left panels 
show the best-fit parameters corresponding to the 
dark matter density profile represented in Fig. [14] 

mates of the PL (nf ) models are significantly larger 
than those coming from the measurements of the 
SIE (nf) models. This is a consequence of the dif- 
ferent error sizes of the two sets of projected total 
mass estimates. For the same reason, as already 
discussed looking at Fig. [Tj|] large values of /l are 
excluded at a 95% CL only if the lens three dimen- 
sional total density profile is fixed to be isother- 
mal. We observe that the dark matter component 
is in any case more diffused than the luminous one. 
In fact, ro /Vl is larger than 2 at more than a 95% 



CL. Given the assumed parametrization, we also 
find that the dark matter density profile pr>(r) is 
probably shallow in the inner galactic regions. The 
value of 7d is indeed lower than 0.7 at a 68% CL. 

5. Summary and conclusions 

By means of HST/ACS and WFPC2 imaging 
and NOT/ALFOSC spectroscopy, we have estab- 
lished that SDSS J1538+5817 is a rare lensing sys- 
tem composed of a luminous elliptical galaxy, lo- 
cated at redshift z\ — 0.143, that acts as a lens 
on two distinct and equally-distant (z s = 0.531) 
sources. The two sources are lensed into a dou- 
ble and a quad (with an almost complete Einstein 
ring) system, covering rather large angular and ra- 
dial scales on the lens plane. This exceptional con- 
figuration has allowed us to investigate in great 
detail the lens total mass distribution within the 
effective radius of the galaxy, through parametric 
and non-parametric point-like lensing programs 
and perform a complete statistical study of the 
errors and correlations on the lens model param- 
eters. Then, by fitting the lens SED with CSP 
models, we have estimated the luminous mass of 
the galaxy and combined the lensing and photo- 
metric measurements to examine the characteris- 
tics of the galaxy dark-matter halo. 

In detail, our main results can be summarized 
in the following points: 

— Parametric models predict image positions 
that match closely the observed lensing ge- 
ometry, and describe lens total mass distri- 
butions that are almost circular in projec- 
tion, moderately steeper than an isothermal 
profile, and well aligned with the lens light 
distribution. 

— The value of the total mass projected within 
the Einstein circle of radius 2.5 kpc is slightly 
larger than 8 x 10 10 Mq and approximately 
10% of this mass is in the form of dark mat- 
ter. 

— In the inner galactic regions, the galaxy 
dark-matter density distribution is shallower 
and more diffuse than the luminous one. The 
former starts exceeding the latter at a dis- 
tance of roughly 6 kpc from the galaxy cen- 
ter, corresponding to 1.5 times the value of 
the luminous effective radius. 
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We conclude by remarking that strong gra- 
vitational lens systems with configurations com- 
parable to or more complex than that of SDSS 
J1538+5817 are excellent laboratories to study the 
distribution of luminous and dark matter in early- 
type galaxies. However, to achieve realistic results 
on the dark matter component, it is essential to 
verify the commonly accepted isothermality of the 
total mass distribution at a higher level than done 
so far. Strong lensing systems with an Einstein 
radius significantly larger than the effective radius 
of the lens galaxy would be invaluable to deter- 
mine the dark matter properties of the halos of 
early-type galaxies. 
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vice mode. We are grateful to the NOT Scien- 
tific Association for awarding some observing time 
solely on the basis of scientific merit and support- 
ing the NOT Summer School where CG gained 
some observational experience with the ALFOSC 
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A. Non-parametric models 

PixeLens ( Saha fc Williams! liool is a non-parametric lensing program that generates an ensemble of 



models consistent with the observed data of a lensing system. Each model is composed of a pixelated surface 
mass density map of the lens, the reconstructed position of the source, and, optionally, an estimate of the 
value of the Hubble parameter. These results are obtained by using the observed positions of the multiple 
images (ordered by arrival time, even if time delays are not known), the redshifts of the lens and the source, 
and some priors based on previou s knowledge of general galaxy mass distribution (for further details, see 
Saha &: Williams] 1997 : Coled 2008 ). Interestingly, PixeLens has been employed to measur e the value of the 
Hubble paramete r from samples of strong lensing systems with measured time delays (e.g., Saha fc Williams] 
2006l : IColesl[2008h . 

We model here the surface mass density of the lens on a symmetric circular grid of 2" radius divided into 
20 pixels. We consider 400 models with fixed cosmological values and with decreasing total projected mass 
profiles [i.e., S(i?) oc R~ a . where a > 0]. 

Total projected mass NP (nf) surface mass density y probability distribution 

- ' i 0.8 1 
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Fig. 16. — Left: Comparison of the total projected mass estimates, with la confidence intervals, from 
parametric [PL (nf)] and non-parametric [NP (nf)] modeling. The arrows show the projected distances 
of the observed multiple images from the lens center. Middle: Isodensity contours of the best-fit non- 
parametric [NP (nf)] total surface mass profile. The observed image positions of the double (diamond) 
and quad (square) systems are shown. Right: Marginal probability distribution (histogram) of the three- 
dimensional total density exponent 7 from non-parametric modeling. The thick bar on the x-axis shows 
the lcr confidence interval. The same probability distribution as obtained by combining strong lensing and 
stellar dynamics measurements is represented by the smooth curves. 



The cumulative total projected mass, the total surface mass density profile of the average model, and the 
marginalized probability distribution of the three-dimensional total density exponent 7 are shown in Fig. 
fTHl We measure a value of the total mass projected within the Einstein radius of 8.59^g \^ x 10 10 M Q , at a 
68% CL. At the same confidence level, we estimate a value of 7 included between 1.62 and 2.87. We observe 
that the contour levels of the non-parametric total surface mass show non-negligible values of ellipticity in 
the inner regions. The differences between the surface brightness of Fig. [TT] and total surface mass of Fig. [T|5] 
are significant within the area defined by the Einstein radius. This is not surprising since here the total 
surface mass density distribution is almost completel y unconstrained by the lensing observables. This is the 
equivalent of Gauss' law in gravitational lensing (see iKochanek 2004 ). We notice that these differences are 
less evident outside the Einstein ring, where the positions of the multiple images limit the freedom of the non- 
parametric models in determining the lens total mass distribution. In Fig. 1161 we also show for comparison 
the mass estimates obtained in the equivalent parametric modeling [PL (nf )] and the probability distribution 
of the density exponent that is expected by combining strong lensing and stellar dynamics measurements. In 
detail, the combined lensing and dynamical probability distribution for 7 is obtained by using the following 
expression 

■,n m ,n A )=(^) 1 s( 7 , <$,/?) (Al) 



An a. 



#Ein 

r(zi,z s 
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that relates through the spherical Jeans equations the values of the central stellar velocity dispersion uq, Ein- 
stein angle #Ein, effective angle 9 e , exponent of the three-dimensional luminosity density profile S, anisotropy 
parameter of the stellar velocity ellipsoid /?, and ratio of angular diameter distances between observer-source 
and lens-source f(zj, z s ] Cl m , Qa) [ff(7 tAi2) * s a numer ical factor that depends on the three cited quantities; 
for definitions and further details, see lKoopmansl[2005| . In the plots of FigfTBl we fix S equal to 2 and choose 
two values of /3 (—0.25 and +0.25) representative of small tangential and radial orbit anisotropy. By doubling 
the size of the grid but keeping the same size of the pixels, we have checked that the choice of a circular grid 
with a radius of 2" to reconstruct the total surface mass density distribution of our not perfectly circular 
lens galaxy does not introduce any artificial shear component and does not affect significantly the results. 

According to these results and looking at Fig. Q1)J we can conclude that the two independent parametric 
and non-parametric analyses are in general consistent, within the errors, as far as total projected mass 
and three-dimensional total density exponent measurements are concerned, but small differences and some 
considerations are worth noticing. 

The projected total mass estimates obtained with PixeLens are systematically larger th an those obtained 
with qravlens. This can be caused by a combination of the mass-sheet degeneracy (see Falco et al. 19851: 
Schneider fc Seitd Il995h and the prior on the positive dcfmitcncss of every pixel of the grid of the total 



surface mass density. Among all the arbitrary constants that can be added to the convergence k, leaving 
though the image positions unchanged, those which provide a negative value of k somewhere on the grid are 
excluded, a priori, from the non-parametric lensing analysis. This fact may bias the projected total mass 
measurements to slightly larger values. 

As far as 7 is concerned, the larger uncertainty coming from the non-parametric reconstruction with 
respect to the parametric one is probably just a consequence of the more general allowed models. A bias 
towards small values of 7 may be associated to the prior present in PixeLens that constrains the value of k 
on one pixel of the grid to be lower than twice the average value of the neighboring pixels. For large values 
of 7, two adjacent pixels located in the central region of the lens may have very different values of k, hence 
these models may not be included in the statistical ensemble. 

Finally, we remark on the overall agreement between the lensing only and lensing plus dynamics probability 
distributions of 7. We notice, though, that lensing alone does not reach the precision needed to distinguish 
among models with different values of the stellar anisotropy parameter (3. 
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